Thermally fluctuating superconductors in two dimensions 
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We describe the different regimes of finite temperature dynamics in the vicinity of a zero tem- 
perature superconductor to insulator quantum phase transition in two dimensions. New results 
are obtained for a low temperature phase-only hydrodynamics, and for the intermediate temper- 
ature quantum-critical region. In the latter case, we obtain a universal relationship between the 
frequency-dependence of the conductivity and the value of the d.c. resistance. 
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In niajW|-|interesting two-dimensional superconducting 
systemstrHa, such as Josephson junction arrays, granu- 
lar superconducting films, and the high temperature su- 
perconductors, it appears that the electrons bind into 
Cooper pairs below a pairing temperature (Tp) that is 
well above the Kosterlitz-Thouless-jxansition, at Tkt, to 
long-range superconducting orderoO'ElLlij. It is natural 
to search for a direction in parameter space where Tkt 
vanishes at a Tj-j^-p superconductor-insulator quantum 
phase transitionO'Ej, and to then expand in the devia- 
tion from the quantum critical point — it is not necessary 
for this point to be experimentally accessible for such an 
approach to be valuable, as it offers a controlled descrip- 
tion of an intermediate coupling regime. We describe 
crossovers in the dynamics near the critical point: new 
results are obtained for a low temperature phase-only 
regime, and for the intermediate temperature 'quantum- 
critical' region. In the latter regime we describe the fre- 
quency (w) dependent conductivity (cr) in terms of a sin- 
gle dimensionless parameter, 7(T'), which determines the 
d.c. conductivity; 7(T) can also be determined by sepa- 
rate static measurements. 

For clarity, we will present our results in the context 
of a familiar microscopic model for the superconductor- 
insulator transition; however, our results generalize to a 
much wider class of systems, and this will be discussed 
towards the end of the paper. We consider a array of su- 
perconducting quantum dots at the sites, i, of a regular 
two-dimensional lattice. The operator hi measures the 
number of electron pairs on dot i, and ipi is its canon- 



ically conjugate phase {[hi 
Hamiltonianliil 



(pj] — Sij). We consider the 



i <ij> 

where A'o is the mean integer number of Cooper pairs on 
each dot, < ij > represents nearest neighbor pairs, Ec 
is the charging energy of a dot, and Ej is the Joseph- 
son tunnel coupling between dots. This model exhibits a 
superconductor to insulator transition as the dimension- 
less parameter g = Ec / Ej is increased through a critical 
value f/c - the phase diagram in the T, g plane is sum- 
marized in Fig |l|. We will be mainly concerned with the 
region g < gc. For T < Tkt, in the superconducting 




9=Ec/Ej 

FIG. 1. Phase diagram of H in (Q). Dashed lines are 
crossover, while the full line is a phase transition. Ps{T) is 
non-zero ior g < Qc, T < Tkt; at T = p^O) ~ {qc - 
with z the dynamic exponent and u the correlation length 
exponent [z = 1 for H). The ground state for g > g^ is an 
insulator. If one studies _ff by a tunneling probe to a good 
superconductor, there will be a zero voltage Josephson cur- 
rent in the region where Ps{T) is non-zero, while for g > gc 
and low T, a Cooper pair current will only flow under a finite 
voltage which overcomes the Coulomb blockade. 



phase, the conductivity has the delta-function contribu- 
tion Re[cr] = {e*^np,{T)/h'^)S{Lj), (e* = 2e), which de- 
fines the superfiuid stiffness Ps{T) (in units of energy). 
The hyperscaling property of the quantum critical point 
implies that all dynamic properties in its vicinity are en- 
tirely characterized by a single energy, Ps(0), which Vjan- 
ishes at g — gc- In particular, the conductivity obeystS 



aiu,T) 



knT' knT 



(2) 



where E is a completely universal scaling function. Such 
a two argument scaling form is not overly constraining, 
and permits a rich variety of behavior which we shall 
describe here as a function of the second argument A = 
Ps(0)/A:bT. Earlier studies of superfiuid dynamics (made 
without reference to quantum phase transitions) emerge 
in limiting regimes, with (|^) only placing restrictions on 
certain parameter values. 

Demanding consistency of (|^) with the defini- 
tion of PstC) above immediately leads to interesting 
conclusionsEj: (i) ps{T)/ ps{0) is a universal function 
of A; (ii) the Kosterlitz Thouless transition occurs at 
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the universal value A = Ac where E is singular; {Hi) 
Tkt = Ps (0) /Ac- Aspects of the data on the high temner|- 
ature superconductors are consistent with such trendsEltl, 
suggesting the proximity of an insulator of Cooper pairs 
(with stripe order). 

We now describe the evolution of dynamic properties 
with decreasing A. 

Phase hydrodynamics: At very large A (fc^T <C Ps(0)) 
vortex excitations are exponentially suppressed, and we 
can derive an effective quantum action, S^p for the contin- 
uum phase variable f{x, t) (where x is a two-dimensional 
spatial co-ordinate and t is imaginary time) in a gradient 
expansion; the resulting action is valid only at the longest 
scales, and is in the spirit of the 'chiral lagrangians' of 
particle physics: 



d^xdr 
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Here /x, are spacetime indices with 9^ = (Vx,9r/w), 
and V is the velocity of the 'spin-wave' excitations pro- 
duced by the harmonic terms in {v remains non- 
singular through the quantum-critical point). The non- 
linear term, arises from integrating out amplitude fluc- 
tuations, and leads to spin- wave scattering; consistency 
of the resulting transport properties with demands 
that the scattering cross-section be universally deter- 
mined by ps(0)-hence the dimensionless number Ai is 
universal. In general spatial dimension, d, Ai multi- 
plies (?it;)2/('^-i)[p3(0)]('^-3)/('^-i), from which it is evi- 
dent that the universality of Ai holds onbi-for 1 < rf < 3; 
we computed Ai for H in an expansior£3 in e — 3 — d 
and obtained 
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Determination of the T > transport properties of 
requires solution of the appropriate quantum kinetic 
equations-we will not do thiS|4iere. The procedure is 
closely analogous to early workllS on phonon transport in 
Galilean-invariant superfluids; an important difference is 
that these systems had a cubic non-linearity which is for- 
bidden in our case by particle-hole symmetry. 
Vortex hydrodynamics: In the shaded region of Fig 
in the vicinity of Tkti the vortices proliferate, and the 
dynaniks can be described by a well-developed classical 
theorylia for the vortices alone. This theory is contained 
within S for A « Ac, and compatibility constrains various 
prefactors. So e.g. for T > Tkt the response of the free 
vortices is controlled by their diffusivity, A2%v'^ / ps{0), 
and their screening length, As^Hv / ps{0))e'^'^/^^''^'^'>'''^ ; 
in general the are arbitrary dimensionless scale 

factors — however, near g = gc they become universal 
numbers. 

Quantum critical: Discussion of this small A regime 
j— l^he remainder of the paper. Previous 



3 — d, or small d — 1, or large N (the number of real 
order parameter components). Here we shall present a 
theory directly for the physical case N = 2, d — 2. 

Unlike the two previous regimes discussed above, it is 
no longer possible to decouple the spin-wave and vortex 
degrees of freedom. We will therefore use the complex 
superconducting order parameter '0(x,i) {t is real time) 
which is the continuum limit of the lattice operator e"^' , 
and allow for both amplitude and phase fluctuations in 
ip. Our theory follows from two hypotheses: (i) The 
equal-time correlations of ip are controlled by a Gaussian 
effective action. Evidence for this hypothesis emerged 
in detailed studies of order parameter correlations in the 
quantum-critical regiont3-the non-Gaussian components 
of the ip correlations are weak because the order param- 
eter anomalous dimension at the quantum-critical point, 
1] w 0.03, is so small. (This small also shows why 
-0 is preferred, over the 'dual order parameter' measur- 
ing vorticitdlj-the latter has an appreciable anomalous 
dimensioned.) [ii] The time evolution of il>{x,t) is de- 
scribed by classical equations of motion. The character- 
istic relaxation (phase coberence) time in the quantum- 
critical regime is of ordeiO h/ksT; however the domi- 
nant spectral weight ip-iS^ frequencies hui < ksT, and 
there is good evidencell3ll3 that the errors made by fo- 
cusing on this low-frequency classical regime are quite 
small. We note that a different classical dynamic model 
was considered in Ref. |l2|, but it does not apply to trans- 
port properties. 

We now define our model for quantum-critical trans- 
port, and then present its exact (numerical) solution. In 
additional to the field 'ip{x,t), we will need the canoni- 
cally conjugate variable measuring density fluctuations, 
Sn{x,t), which is the continuum limit of fii — Nq. The 
equal-time correlations of ip and 6n are described by the 
partition function 

Z = J VtP{x)VSn{x) exp(-(Hi + H2)/A:bT), (3) 

Hi= I d'x[\Wi^\'+rHTM^], (4) 

n2 = (l/[2xu(T)]) / dMSn{x))'. (5) 
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relied upon expansions in either small 



The overall scale of ip is arbitrary, and has been adjusted 
to obtain a unit coefficient for the gradient term in (|^). 
The correlation length, f(T), (assumed to be larger than 
all microscopic scales), and the compressibility Xu{T), 
are determined by the underlying quantum physics, and 
will generally be regarded as unknown input parameters 
in our dynamic theory; proximity to the quantum-critical 
point does (weakly) restrict their allowed T dependen- 
cies, but we can evisage a more general application of 
our dynamic theory, without reference to quantum criti- 
cality, in which case ^(T), XuiT) are arbitrary. The time 
evolution of 0, 6n, is defined by a minimal model — the 
Hamilton- Jacobi equations of Tii +Ti.2, and the only non- 
zero Poisson bracket 
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{6n{x),i:{x')} 



P.B. 



{i/n)'il;{x)S^{x-x'), 
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which is the continuum, classical limit of the commu- 
tatator between ipi and fii. Notice that h appears on the 
r.h.s. even though we are considering classical equations. 
The equations of motion implied by (Q)-(|^) are simple 
and familiar. They are the continuity equation 



dSn{x,t)/dt + S/ ■ J{x,t) = 0, 
where J = — (z/?i)('0*V'0 — ipVij)*), and 

dTp{x,t)/dt = -[i/K)^{x,t)^{x,t), 



(7) 



(8) 



which is the Josephson equation, with the electrochem- 
ical potential ^{x,t) — 5n[x,t)/xu{T)- We restate our 
dynamical model for quantum critical transport: choose 
a set of equal-time initial conditions for and 5n from 
the thermal ensemble defined by Z, and evolve them de- 
terministically by and (||). Correlation functions are 
determined by the average over initial conditions. 

We will shortly present convincing numerical evidence 
that define a sensible continuum theory free of 

short distance or short time ('ultraviolet') divergences; 
this is supported by perturbative and renormalization 
group arguments, which we do not describe here. So un- 
like Sip, the present theory describes the couplings of vor- 
tex and spin-wave fluctuations at different length scales. 
It is helpful to visualize the continuum theory by coarse- 
graining to a lattice spacing, a: the shorter distance de- 
grees of freedom lead to fluctuations in the phase and am- 
plitude of ■0, but the long distance transport properties 
are insensitive to the value of a. The absence of factors 
of a in the final results allows us to deduce their func- 
tional dependence on ^(T), Xu{T) by simple engineering 
dimensional analysis. In this manner we conclude that 



a(a;,r) = — 7(r)S, 



(9) 



where Sc is a universal function ('exact' numerical results 
are below), and the dimensionless ^{T) is defined by 



7(r) ^ [kBTxu{T)e{T)f'^. 



(10) 



Consistency of (^) with ^ only requires that 7(T) is 
a universal function of A (also recall that the present 
quantum-critical theory is valid only for small A). The 
result (^) has a clear experimental signature: it implies 
a correlation between the value of the d.c. conductiv- 
ity and the inverse frequency-width of cr(ti'), as they are 
both determined by j(T). More stringent comparisons 
can be made by using the measured d.c. conductivity to 
determine the unknown 7(T), and then using our numer- 
ical results below for Ec to determine the lo dependence 
of (7. It should be noted that our present computations 
for Sc will not apply for very large hLo/ksT, for then a 
full quantum theory is necessary, and wfl||e*ossover to the 
phase-coherent regime discussed earliercHla. 



Exactly at g = Qc, hyperscaling arguments imply that 
£_{T) ~ T-i/^ and Xu{T) ~ T^^-^)/^; so j{T) is T- 
independent for any z, and is expected to be a universal 
number. For thejaodel in (^, this universal number 
can be computedlla in a 1 /N expansion: 
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We turn to our numerical results for Ec. The sim- 
ulations were carried out on an A'^ x square lattice 
of spacing a with periodic boundary conditions, and 
a ^ £,{T) ^ Na. For each a, successively larger 
values of N were used, until the results became iV- 
independent, and the continuum limit was then ap- 
proached as a — > 0. The lattice form of Hi + 7^2 was ob- 
tained by mapping the momentum (fc^,, ky) dependence of 
the couplings of the continuum Hamiltonian under k^ — > 
{4:/a'^)[Ki sm^{k^a/2) + K2 sm^{k^a) + K3 sm'^{Sk^a/2)], 
and similarly for ky (for |V'(/'P, this amounts to including 
first, second, and third neighbor couplings between the 
^p). We chose two different sets of values for if 1^2, 3 to test 
the independence of the continuum theory on the lattice 
realization — (A) the familiar Ki = 1, K2 = K3 = 0, 
and (B) Ki = 3/2, K2 = -3/20, K3 = 1/90 for which 
kl fc^(l -I- O(fc^a)^). The initial conditions specified 
by Z form a Gaussian ensemble, and are easily generated 
in a single sweep. The time evolution was carried out by 
a fourth order predictor-corrector algorithm with a time 
step determined to conserve total energy to a relative 
accuracy better than 10~^. We measured the autocor- 
relation function of the total current, Cj, by averaging 
over 3000 initial conditions; Cj is normalized such that 



Ee(cj)= / dtcj{t)e 

Jo 
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where uJ = j{T)nuj/ {ksT) and t = keTt/ Ck-yiT)). 

Our results for Cj and Sc are in Fig Notice the 
a dependence at t = — this is expected as an ana- 
lytic calculation of equal-time correlations shows that 
Cj(0) = (I/tt) ln(l/a) + . . .. However, the a dependence 
quickly disappears at small non-zero t, and then a uni- 
versal continuum value obtains; this is strong evidence 
for the existence of the continuum theory (p|)-(H). We 
obtained Ec(0) « 1.85; for the model H in (|l|j^ this com- 
bines with (|l|) to yield cr(0,T) = 0.82LE*'^/h) at g = g^. 

The experiments of Rimberg et ala appear to be a 
convenient testing ground for our theory. The metallic 
gate scijeens the long-range Coulomb interactions and so 
justifie£3 the short-range coupling in 7^2- However, it 
also introduces dissipation, which will induce additional 
terms in the equations of motion e.g. (H) is modified to 



d^/dt = -T5ni/6i;* - {i/h)<S>'ilj + c. 



(13) 



where F is a damping co-efficient and ( is the associ- 
ated random noise. Fortunately, it can be shown, as in 
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FIG. 2. Current autocorrelation, Cj, of (|3[)-(|8|), with lattice 
realization A, iV = 64, and a/^{T) = 1/8 (crosses), ^2/16 
(open circles), 1/16 (triangles). Inset: real and imaginary 
parts of Ec, obtained from (lid) for a/^(r) = 1/16; for large 
|t<J| we have ReEc ~ ^/\^\ and ImEc ~ ln(|£l;|)/Zj. 



Refs g^, that for large ^(T), such modifications leave the 
dynamic function Ec unchanged. The damping can be 
a releiiant perturbation on the T = quantum critical 
pointtil, but this manifested only via changes in £,{T), 
Xu{T), and 7(T). (Outside the quantum critical region, 
the change in quantum universality leaves the form (|^) 
unchanged, but with a new E; the structure of the phase 
hydrodynamics, S*^ can also be modified. Similar com- 
ments apply to the following paragraph.) 

As a second modification, we consider the inclusion 
of long-range Coulomb interactions. This will modify 
7^2 to e*2/ tPx(fx'Sn{x)6n{x')/2\x-x'\, and leave ^) 
unchanged, but with<i>(x,t) = e*^ / d'^x' Sn{x' ,t)/\x—x'\. 
Assuming the existence of the continuum limit of the 
classical dynamic theory, we again obtain (^), but with 
7(r) = [fcBT^(r)/e*2]i/2; the functional form of wih 
of course be changed and requires a separate numerical 
simulation. Note that at g = gc, ^{T) is a T-independent 
universal number only .-if z = 1, the value expected for 
long-range interactionsH. j-. 

The experiments of Corson et alu on an underdoped 
high temperature superconductor observe scaling closely 
related to (||), with the frequency scale proportional to 
the resistivity scale, and a scaling function with large 
Lo behavior similar to that in Fig y. However their fits 
use a prefactor Tg whose weak T-dependence disagrees 
with our continuum two-dimensional theory. We spec- 
ulate that this can be explained by corrections to scal- 
ing which are appreciable because both modifications dis- 
cussed above are present here - damping from fermionic 
excitations and long-range Coulomb interactions. Corson 
et aZ.cl motivated the scaling using the vortex theoryliJ, 
but did not consider bound vortex pairs, whose contri- 
butions do not obey their scaling assumptions. 

In summary, this paper has delineated the distinct 
dynamic regimes of thermally fluctuating superconduc- 
tors: the low T phase-only hydrodynamics, the classi- 
cal vortex hydrodynamics in the vicinity of Tkt, and 



the quantum-critical region where phase and vortex fluc- 
tuations strongly coupled — here we proposed a dynami- 
cal model in which non-linear 'mode-coupling' terms de- 
manded by the Poisson bracket (|^) dominate the univer- 
sal, low- frequency, dissipative dynamics. Our approach 
could be extended to other two-dimensional systems: (i) 
the magnetic field-tuned superconductor-insulator tran- 
sition, where a coupling to the external field would be 
required in (||), and [ii) quantum Hall transitions. 
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